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Abstract 

This paper shows how mesoscopic nonequilibrium thermodynamics can be ap- 
plied to condensation and evaporation. By extending the normal set of thermody- 
namic variables with two internal variables, we are able to give a new theoretical 
foundation for a mechanism of condensation that has been proposed from molecular 
simulation results. The flux does not follow a simple Arrhenius formula for small 
activation energies which are relevant here. 

1 Introduction 

Evaporation and condensation has been well studied, and many theories exist for this 
important phenomenon. The theories are mostly rate theories. Kinetic theory has long 
played a central role £Q. With information about the condensation coefficient in this 
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theory, the mass and heat fluxes are determined. Recently Nagayama and Tsuruta j2j 
have given a derivation of the condensation coefficient using transition state theory 
They based their theory on molecular dynamics simulations of condensing particles jlj, 
giving that the probability of a particle to condense increases with an increasing normal 
component of the velocity of the particle to the interface. 

Common to rate theories [Sj is that the driving forces of transport do not appear 
explicitly in the expression for the rate. Nonequilibrium thermodynamics gives, on the 
other hand, two independent forces for the case of condensation or evaporation of a one- 
component system. These are the chemical potential difference at constant temperature, 
and the difference in inverse temperature jSlEl- I n the nonequilibrium thermodynamic 
description of phase transition, the heat and mass fluxes are linear in these forces. 

Rubi and coworkers |H| developed mesoscopic nonequilibrium thermodynamics, fol- 
lowing and generalizing the work described in de Groot and Mazur This method, 
which describes the system on the mesoscopic level, can give a more detailed thermo- 
dynamic foundation for rate processes. This is done by including more variables than 
are normal in a nonequilibrium thermodynamic description, for instance the velocity of 
the particles, an important variable in rate theories. By integrating the flux equations, 
which are now linear in the forces on the mesoscopic level, nonlinear flux-force relations 
are obtained. 

Rate theories and nonequilibrium thermodynamics aim to describe the same phe- 
nomenon. It is therefore important to know the strength and weaknesses of the theories, 
and under which conditions they become equivalent. It is the main aim of this article to 
develop and apply mesocopic nonequilibrium thermodynamics to evaporation and conden- 
sation. By developing the theory for this case, we hope to strengthen the thermodynamic 
basis for this important phenomenon. We shall present a detailed mesoscopic formulation, 
starting with the macroscopic one, in order to define proper boundary conditions. As new 
variables we introduce the condensation coordinate, called the reaction coordinate by Na- 
gayama and Tsuruta, and the normal velocity of the particles as internal variables. By 
doing this we shall see that the earlier description using nonequilibrium thermodynamics 
for interfaces El can be related to a richer, more detailed formulation. The new level 
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of description is also easier to compare to kinetic theory. This theory will appear then as 
a special case. 

The paper is organized as follows. The nonequilibrium theory for the mesoscopic level 
is developed in section 2. We shall use the terminology used in de Groot and Mazur jU]. 
Boundary conditions are given for the fluxes along the internal coordinates. The result of 
the section is a new expression for the excess entropy production rate. In the 3rd section 
the description is restricted to stationary evaporation and condensation. This simplifies 
the expression for the excess entropy production rate considerably. The resulting flux 
equations are next given in section 4. For large enough times one may neglect inertial 
effects due to acceleration of particles. This leads to a simplified expression for the flux 
along the condensation coordinate. Leaving the effects of the thermal force to future 
work, we then proceed to neglect temperature differences across the interface. Integrating 
the resulting expression for the flux along the condensation coordinate over the internal 
variables, an expression is found for the condensation flux across the interface in the 
isothermal case. Our result for the condensation flux is discussed and compared to the 
corresponding ones found using simple kinetic theory and in the work by Nagayama and 
Tsuruta [2], in section 5. They are found to reproduce the results of both these cases. In 
the final section some conclusions are drawn. 



2 The system on the macroscopic level 

The system consists of a vapour in the presence of its liquid in a container of volume 
V. The interface is planar. On a macroscopic scale (say in /iin), the vapor is located 
in the region with x < 0, while the liquid is located in the region with x > 0, see 
Fig. la. The interface appears as a discontinuity on this scale, but with excesses in entropy, 
internal energy, and mass. This formulation is due to Gibbs PU] • I n the absence of global 
equilibrium, there are gradients in temperature, T, and density, p, in the homogeneous 
liquid and vapour phases. Across the interface, there is also a local variation in the same 
quantities, as well as in the velocity distributions of the particles. This variation is, as 
Nagayama and Tsuruta [2] have shown, important for condensation. 
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The expression for the excess entropy production rate of the interface, in the presence 
of a temperature difference and a difference in chemical potential was given by Bedeaux 
and Kjelstrup [TT]. For a one-component system we have: 
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The flux of heat from the vapour into the interface is J'' 9 , while the one from the interface 
into the liquid is J 1 ' . The corresponding mass fluxes are J 9 and J 1 . When the independent 
fluxes are the measurable heat fluxes and the mass fluxes on the vapor and the liquid sides, 
the conjugate forces are the differences in inverse temperatures (from the interface to the 
vapour, and from the liquid to the interface) and the same differences for the chemical 
potentials at the temperature of the interface divided by the temperature of the interface. 
The superscripts g, s and I indicate values in the vapor, at the interface and in the liquid, 
respectively. The chemical potential in the interface is /i s . When we write fi 9 (T s ) it 
means that the chemical potential of the gas is taken at the temperature of the interface. 
When the chemical potential of the gas is taken at the temperature of the gas this is 
not separately indicated, thus n 9 = /i 9 (T 9 ). The frame of reference for the fluxes is the 
interface. For simplicity we do not consider fluxes along the interface in this paper. 



3 Mesoscopic nonequilibrium thermodynamics 

The purpose of this work is to provide a mesoscopic description of the process, in terms 
of internal variables characterising the transition states of the system, able to derive the 
kinetic law for the rate of condensation. 

We shall replace the third and the fourth terms in Eq. (Q) by a more detailed expression 
derived using two internal variables. The choice for internal variables is motivated by 
the work of Nagayama and Tsuruta [2|. We know from their work that the potential 
energy of condensing particles has a maximum in the interfacial region and that the 
velocity component normal to the interface is important for the likelihood of condensation. 
It is also known that the energy flux across the interface has maximum contribution 
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from potential energy in the interface [12]. In order to describe this, we introduce a 
coordinate 7 for the degree of condensation. This coordinate is similar to the degree of 
reaction in transition state theory for chemical reactions [3]. By introducing the degree of 
condensation as a variable, one enlarges the space of the interface into an internal space, 
ranging from 7 = to 7 = 1. This is illustrated in Fig. lb. In other words, we expand 
the interface in this coordinate (at x = 0). In the process, we must make sure that the 
boundary conditions are correct. 

In FigQk the chemical potential has a separate value in the interface. The problem 
is to find the chemical potential in the interface as a continuous function of the internal 
coordinate, i.e. to find the function /i s (7). We need boundary conditions for the chemical 
potential that are compatible with Eq.([Tj). The boundaries for the chemical potential 
^(7) in the 7-space are naturally: 

»' = m wd m = £ (2) 

Additionally, the flux along the condensation coordinate, J* (7), must obey 

j;(0) = J 9 and J s 1 (l) = J l (3) 

We shall also use as internal coordinate, the component of the velocity normal to the 
interface, v which varies from v = —00 to v = 00. The chemical potential will then be a 
function of both 7 and v. The relation between ^(j) and /i s (7,v) will be given in the 
next section. The flux along the internal coordinate v, J*, is more commonly known as 
the acceleration of the particles. It satisfies the following boundary condition: 

J*( 7 ,-oc) = and J v *( 7 ,oo) = (4) 

We furthermore have 

/oo 
*i;(7,v) (5) 
-00 

This means that the value of J* integrated over all velocities of the condensation coordi- 
nate is equal to the mass flux from the vapor into the interface. There is no contribution 
to the mass flux at infinitely high or low velocities at any coordinate 7. 
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3.1 The chemical potential in internal variable space 

The variation in the excess entropy of the interface, s s , in this mesoscopic description is 
given by 

Ss s = ^6u s - -L J d 1 J°° dv //( 7 , v)tf/'( 7 , v) (6) 

The excess entropy of the interface is a function of the excess internal energy of the 
interface, u s , and the probability, / s (7,v), to find a mole of particles in the interfacial 
area with the internal variables 7 and v, see [S|. As the equation shows, the derivatives 
of the entropy with respect to these variables are the inverse interface temperature, 1/T S , 
and the chemical potential, /z s (7,v), divided by the interface temperature, T s . In the 
general case, with / s (7, v), there is a probability for the particle to have a particular 
velocity The velocity dependence is integrated out: 

/CO 
dv/'(7,v) (7) 
-co 

to give the probability distribution / s (7). This distribution can be interpreted as the 
probability to find a mole of particles in the interfacial layer in Fig. lb with the internal 
variable 7 (in mol/m 3 ). 

According to Gibbs entropy postulate P, we can write that 

/ i s (7,v)= J RT s lnr( 7 ,v) + C s (7,v) (8) 

The potential function C s (7,v), which depends on 7 and v, may also depend on T s . 
In equilibrium the entropy has a maximum for a given value of the temperature and it 
follows from Eq.lJHjl. that the chemical potential is constant. We shall use the constant 

/4(7, v) = H*(T eq ) - RT eq In (^^) (9) 

The first term to the right, /j,*(T eq ), is the chemical potential of the saturated vapor, 
considered as the standard state. This chemical potential is independent of 7 and v. The 
last term is minus two times the translational potential energy of one mole ideal gas. 

Alternatively, we can use the Gibbs postulate for the description where the velocity 
coordinate has been integrated out. We then have 

/ i s ( 7 ) = J RT s lnr( 7 ) + ^(7) (10) 
6 



In equilibrium, we have 



The second contribution is half the value of that in Eg. ([§]). The relation between the two 
expressions for the chemical potentials, Eqs.(J0J) and is found after determining the 
values of the functions C s (7, v) and C s (7), see Eq. fjTTj) below. Those values C s (7, v) and 
C s (7) will be obtained from the probability distributions at equilibrium between liquid 
and vapour. 

Consider first the determination of C s (7, v). The probability varies with the kinetic 
energy of the particle which for a mole of particles, is given by Mv 2 /2. We shall further 
assume that the particles are crossing a potential energy barrier along the condensation 
coordinate, 0(7). The barrier is zero at the ends of the condensation coordinate, 0(0) = 
0(1) = 0. In equilibrium, the probability distribution is then given by: 



A(7,v) = exp 
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By substituting this equation and Eq.© into Eq.(JHJ), we find 



C s (7, v) = 0( 7 ) + -My 2 



(13) 



The function C s (7, v) is illustrated in FigJ21 This part of the chemical potential has 
a small maximum along the 7-axis. Particles at the left of the barrier having higher 
velocities (and consequently higher values of C s ) can easily cross the barrier in 7-space 
going through decreasing values of the velocity coordinate. This shows that the velocity 
distribution of the particles plays an important role in condensation kinetics [2]. 
On the other hand, the probability distribution in the 7-space is: 

>eg(7) -0(7)' 



A (7) = ex P 
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By substituting this equation and Eq. ffTTj) into Eq. ffTUJ) . we find 

C s (7) = 0(7) 



(14) 



(15) 



The expressions are now in accordance with each other. This can be seen by first 
substituting / s (7,v) and / s (7), as obtained from Eqs.(jHI) an d (jTOjl . into Eq.(J7J). The 



following identity is obtained: 
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Dividing both sides by exp[— C s (7) / RT S ] , and using Eqs.([13|) and f|15j) . this can be 
written as 



cxp 



^(7) 



dv exp 



y( 7 , v )- iMv 2 



(17) 



Substitution of Eqs.(|9j) and (jTTj) into this equation gives an identity This relation explains 
the factor 1/2 in front of the last term in EqJll| while there is a 1 in the equation before 
it. 

The boundary values for the chemical potential in 7-space, can now be found using 
Eq. pij) and the expression for \x s . We obtain: 
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The first equation says that the states with 7 = and 1 are identical at equilibrium, 
whereas the second equation assumes ideality. 

In summary, we have given above two sets of chemical potential and corresponding 
probability distribution. In the most general case, the chemical potential is a continuous 
function of 7, v. In the other case it is a function of 7. The relationship between /x( 7 , v) 
and /x( 7 ) was given by integrating the most general probability distribution function 
over v-space. An important consequence of the change in variables, is that the chemical 
potential in the state of equilibrium changes. As the number of internal variables increases 
the equilibrium state changes to lower values. With one internal variable, the chemical 
potential is lowered by half the kinetic energy of an ideal gas, with two internal variables 
it is lowered by twice this value. 
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3.2 The entropy production rate as a function of internal vari- 
ables 

We can now follow the normal procedure and find the entropy production rate in the 
variable space we are interested in. The procedure consists of combining the time rate of 
change of Gibbs equation with the balance equations. In the present case, there are three 
independent balance equations. 

The balance equation for the probability distribution is 

^ = -^-)-^,v, (19, 

By integrating this equation over the velocity coordinate, we obtain the balance equation 
for probability distribution along the condensation coordinate alone 

The conservation equation for internal energy is 

%■ = J? - A (21) 

where each energy flux, Jf or Jl is a sum of the measurable heat flux and the co-moving 
enthalpies, J' q + hJ, on the vapor side and on the liquid side of the interface, respectively. 
The balance for the entropy of the interface is 

%-V-m-o- (22) 

where the entropy flux is given by J s = (J e — /iJ)/T = J'/T + sJ in both bulk phases 
(with superscripts g and /, respectively). 

The time rate of change of the interfacial entropy density is obtained from Eq.©: 

dt T s dt T s J ' J-oo & 
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By substituting Eqs. (fT9|) and (j2Tj) into this equation, we find 
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^7 / dv J*( 7 ,v) — /i s (7,v) 
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(24) 



where we have used partial integration and the fact that J* is zero for v = ±oo. Comparing 
this equation with Eq. (J22|) we may identify the entropy production rate 
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In order to simplify this equation we need boundary conditions for 7 = and 1. For 
7 = the particles leave the vapor and start moving along the condensation coordinate. 
At 7 = 1 the particles are at the end of the condensation coordinate and enter the liquid. 
In the vapor and the liquid we do not use the velocity coordinate as an internal variable. 
This corresponds to the assumption that the chemical potential in those regions does not 
depend on the velocity. At the ends of the condensation coordinate, we shall similarly 
assume that the chemical potential is independent of the velocity. In agreement with our 
previous choice for standard states, we have using Eq. (fT7J) : 



fi s (0,v) + ^RT s \n 



( 2nRT t 
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//(l) (26) 
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Substitution of Eg. (1261) into Eg. (1251) results in 
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The first two terms are identical to the first two terms in Eg. ((J). When /i s ( 7 ,v) = /j, 8 is 
constant, the last two terms are zero and the fourth and fifth term are identical to the 
third and fourth term in Eg.([Tj). The third term is not independent of the two following. 
It is zero in the stationary state, however. The last two terms provide the new element 
in this mesoscopic description. 



4 Stationary evaporation and condensation 

Most experiments, Monte Carlo and molecular dynamics simulations are performed for 
stationary evaporation and condensation. As the description of this case simplifies con- 
siderably, we shall restrict ourselves to this case. For the stationary case it follows from 
the balance eguations that 



(28) 



(29) 



Jf = J l e = J e and J 9 = J 1 = J 
For the measurable heat fluxes, this relation gives 

J'/ = J' q 9 + J(h l - h 9 ) 
Furthermore it follows from Eq. (jl9|) . using stationarity, that 

|^(7,v) + |^( 7 ,v)=0 (30) 
Integrating this equation over v, using J® (7, —00) = J* (7, 00) = and Eg.(JSJ), gives 

rfvj;( 7l y) = j;( 7 ) = J (31) 
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Using these relations the excess entropy production rate may be written as 
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where we used Eqs. (j2*U|) and (}3*T]) . In the last identity we have also employed the thermo- 
dynamic identity 
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where the differentials are evaluated at constant pressure. The pressure is constant be- 
cause the system is in mechanical equilibrium. 

If one uses the energy flux rather than the heat flux in Eq. (j32)) it becomes 
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In this form it is clear that in order to transfer a particle from the vapor to the liquid, it 
must first be transferred from the vapor to the beginning of the condensation coordinate, 
then it has to go along this coordinate to its end and finally the particle should be 
transferred to the liquid. The total resistance to go from one end to the other is the sum 
of the resistances for these three steps. The description without internal coordinates is 
obtained by setting /I s (7, v) = fi s in Eq. (p£fj) . In that case only the first two terms remain. 
The force giving mass transfer has simplified to fi l /T l — [i 9 /T 9 . We see that the meaning 
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of introduction of internal variables is to distribute this force along the condensation 
coordinate. In line with the introduction of the internal coordinates Eq.(j2J, is a natural 
choice to make. This then simplifies Eq. (j34j) to 
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5 Flux equations 



The stationary state excess entropy production rate that was derived above has three 
fluxes and forces. Nonequilibrium thermodynamic theory prescribes that each flux is a 
linear combination of each force of the system. Using Eq. (J34|) . the linear relation for the 
fluxes in terms of the forces become 
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J e = L ee A(—)- I c?7(ivL e7 (7,v) — — ^ / d 7 dvL ev ( 7 , v) 
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J 7 ( 7 ,v) = L 7e (7,v)A(-) -L 77 ( 7 ,v) — — ^ L 7V ( 7 , V ) — ^ 

J v (7,v)=L ve ( 7 ,v)A(i)-L v7 ( 7 ^ (36) 

where A(l/T) = 1/T — 1/T 9 . All conductivity coefficients except L ee are functions of 
the internal coordinates, while they all depend on T s and c s . The matrix of conductivity 
coefficients satisfies the Onsager symmetry relations. 



5.1 Mechanical equilibrium 

A non-zero flux of particles in velocity-space, means that there is an acceleration, which 
is partially balanced by inertial effects. In the absence of acceleration, J v =0, and the 
system is in mechanical equilibrium. Equation (J55)l gives: 

M7,v)A(I) _ ivl(7 , v) ^^Z) _ ivv(7 , v) |.^) = o (37) 

The last term, that contains a main coefficient, is a dissipative term (the intertial term), 
while the two first terms represent work terms. The central term is the work done by 
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moving a mole of particles up the barrier in the internal variable space. We can use this 
relation to eliminate dfi s ( , y,v)/dv from the other fluxes. This gives for the flux along the 
condensation coordinate 



J 7 ( 7 ,v) = L 7e ( 7 ,v)A(i) - L 77 ( 7 ,v)|-^f^ (3? 



where 
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5.2 The isothermal case 

When the system in mechanical equilibrium is also isothermal, Eq.(|38|) further reduces to 

J 7 (7,v) = -L 77 (7,v)^g^ (40) 

where T is the uniform temperature. It is now convenient to write the gradient of the 
chemical potential in an alternative form. Using Eqs.(jHJ) and (|13|) . we have 

//( 7 ,v) = iff* In ( 7) v) + 0(7) + \mv 2 = iff In [f (7, v)/ 7 s ( 7 )/ v s (v)] (41) 
where we defined 

it!Tln/ 7 s (7) = 0( 7 ) and iff In /» = l -Mv 2 (42) 
Using now Eq. (|4*TJl we write Eq.(|40|) as 

7(7 ' } " Hi^mmd!^ rt (43) 

It is commonly found that the conductance is proportional to / s (7,v), so we will assume 
that the diffusion coefficient along the condensation coordinate, which we define by 

/ s (7,v) 

is independent of 7. Equation (j43|) can then be written in the form 
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Integration of this equation over 7 gives, using Eq.(j2j), 
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Further integration of this equation over v, taking Eq. (J!?Tj) into account, yields: 
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For the chemical potential of the vapor we have, using Eqs.(j2j) and ( I18)). 
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The chemical potential of the liquid corresponds to the chemical potential of the saturated 
vapor at the temperature of the liquid: 



H l = ^(1) = fj,* (T) - -RT 'In 
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Substitution of these expressions in Eq. (j4Tjl . and using Eq.(|18j) results in 
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We define the first integral to the right as 



exp 



RT 



1 -1-1 



LJQ 



(51) 



where we E Q is an activation energy. The final result for isothermal phase transition at 
mechanical equilibrium is then: 
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The flux of mass according to this equation is driven by the difference between the 
vapour pressure and the saturation pressure. It is proportional to an integral that depends 
exponentially upon the velocity of the particles that are hitting the interface, as well as a 
diffusion coefficient in velocity space. The flux depends also on an activation energy i?o, 
in a common way. 
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6 Discussion 



The theory we have presented can in general be applied to cases that are away from 
mechanical equilibrium and are not isothermal. Under mechanical and isothermal condi- 
tions it gives expressions for the condensation rate which agree with those obtained from 
kinetic theory and transition state theory. We will proceed in this section to show that 
agreement and will leave for the next section a general discussion to summarise the main 
characteristics and to point out some possible extensions. 

6.1 The simplest version of kinetic theory 

Kinetic theory provides expressions for the unidirectional rates of evaporation or conden- 
sation, and the net rate as a difference between these. It establishes that the evaporation 
rate is proportional to the pressure of the saturated vapour (which gives the state of the 
liquid), while the condensation rate is proportional to the actual vapour pressure PQ. This 
gives as maximum condensation rate for isothermal conditions: 

J = (2vr MRT)~ 1/2 [p - p*(T)} (53) 

In the most simple kinetic theory, there is no potential barrier along the reaction 
coordinate (E = 0). In our formulation this means that /^(j) = 1 in our expressions. 
We then reproduce the result from kinetic theory by using: 

D(v) = v (54) 

Substituting this into Eq. (J3UJ) and restricting the integration to positive velocities one 
obtains Eq.([S3J). 

If one takes a potential barrier into account, as we will do below, we obtain a prefactor 
in Eq. (p)3*j) . The prefactor, which is the condensation coefficient, gives the fraction of 
incident particles that condenses. The equation then becomes equal to the Hertz-Knudsen 
formula for isothermal conditions [T3JE]- in Schrage's formulation jTJ], the Maxwellian 
distribution of the incident particles is corrected by a possible net bulk velocity. This type 
of correction will not be considered now. Neither will we consider more recent results of 
kinetic theory. 
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6.2 Modified transition state theory 

Nagayama and Tsuruta [2] found that the condensation coefficient depended on the in- 
cident velocity of the particles. Their system was in equilibrium and had no net bulk 
velocity. They calculated the activation energy as the difference in the mean kinetic 
energy of condensing and Maxwellian particles. 

If there is a potential barrier along the condensation coordinate, we have ^(7) 7^ 1. 
This also affects the expression for the diffusion coefficient D(y). The expression for the 
condensation flux then becomes 



J = a c (T)(27rMRT)- 1/2 \p-p*(T)\ 



(55) 



where the prefactor, the condensation coefficient, is given by 



a c (T) = exp 



E, 



, MD(v) 

, , dv exp 

RT J J RT 1 



Mv 2 ' 
2RT 



(56) 



The dependence of the fraction of condensing particles on the incident velocity j^j, 
modifies the expression for D(v) to 

Mv 2 ' 



D(v) = v l-/3exp 



For the condensation coefficient we then obtain 



2RT 



(57) 



(5f 



Following Nagayama and Tsuruta |2j, the activation energy is equal to the average kinetic 
energy of the condensing particles minus the average energy they would have for (3 = 0. 
This gives 



En 



f dv \Mv 2 v (1 -(3 exp 



Mv 2 
' 2RT 



exp 



Mv 2 
2RT 



J-dvv (l-^expI-fg^expHg] 



Jo°° ^ v ^Mv 2 v exp 



Mv 2 
2RT 



f °°dvv exp 
Substitution of this equation in Eq.(|3S|l results in 



1 P 1 
1 - - J exp 



RT 







2(2-/3) 



(59) 



2(2-/3) 



(60) 



17 



which correponds to the result Nagayama and Tsuruta j2] found. These authors go on to 
give a further expression for (3 in terms of the system's specific volume ratios. One can 
also express the condensation coefficient in terms of the activation energy: 



The flux does thus not follow a simple Arrhenius formula, not even for small E . The 
deviation from Arrhenius behaviour is illustrated in FigOl The figure shows that the 
value of the condensation coefficient is always significantly smaller than expected from an 
Arrhenious formula. Nagayama and Tsuruta found Eo/RT < 1. The importance of these 
expression is that they make it possible to predict a priori the condensation coefficient 
(or the flux) from the system's equilibrium properties. 

7 Conclusions 

In this paper we have proposed a new theory to deal with evaporation and condensation 
phenomena in nonequilibrium situations. Our theoretical framework is based on meso- 
scopic nonequilibrium thermodynamics 0|B|, a theory proposed to explore the mesoscopic 
level of description, using the tenets of nonequilibrium thermodynamics. 

We have derived an expression for the condensation flux, which can be applied to non- 
isothermal situations or even to the case in which the system is away from mechanical 
equilibrium, when acceleration of particles becomes significant. 

Our result reproduces the ones obtained using different approaches. In the absence 
of a potential barrier along the condensation coordinate, we rederive the expression for 
the rate given by the simplest version of kinetic theory. Using the expression for the 
condensation probability obtained by molecular dynamics of condensation, we obtain the 
condensation coefficient found by Nagayama and Tsuruta 

Our theory has shown how the method, set up by nonequilibrium thermodynamics |9] 
to deal with transport processes, also holds for activated processes. In an early work of 
Eyring and Eyring [3| that possibility was questioned. An unified formalism, able to cope 
with the difficulties inherent to nonequilibrium situations, has then been established. 




(61) 
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In this sense, we can elucidate nonequilibrium situations in which the system as a whole 
is subject to temperature and/or velocity gradients. These are interesting generalisations 
of earlier conditions used. We hope in this way to contribute to the understanding of 
phenomena which pose interesting theoretical problems and are practically important. 
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Figure Legends 



1. The interface separating the homogeneous gas and liquid phases on a /im-scale (a) 
and in internal variable space (b). 

2. The potential function C s (7, v) as a function of the internal variables 7 and v = 
v/(2RTM- 1 )- 1 / 2 . 

3. The condensation coefficient as a function of the activation energy (solid line). The 
dotted line corresponds to Arrhenius behaviour. 
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